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Mechanical shear deformations lead, in some cases, to effects similar to those resulting from ion 
irradiation. Here we characterize the effects of shear velocity and temperature on amorphous silicon 
(a-Si) modelled using classical molecular dynamics simulations based on the empirical Environ- 
ment Dependent Inter-atomic Potential (EDIP). With increasing shear velocity at low temperature, 
we find a systematic increase in the internal strain leading to the rapid appearance of structural 
defects (5-fold coordinated atoms). The impacts of externally applied strain can be almost fully 
compensated by increasing the temperature, allowing the system to respond more rapidly to the 
deformation. In particular, we find opposite power-law relations between the temperature and the 
shear velocity and the deformation energy. The spatial distribution of defects is also found to 
strongly depend on temperature and strain velocity. For low temperature or high shear velocity, 
defects are concentrated in a few atomic layers near the center of the cell while, with increasing 
temperature or decreasing shear velocity, they spread slowly throughout the full simulation cell. 
This complex behavior can be related to the structure of the energy landscape and the existence of 
a continuous energy-barrier distribution. 



I. INTRODUCTION 

The structural and dynamical properties of amorphous 
silicon (a-Si) have been extensively investigated using 
classical fo rce field s^^, tight-binding^^ and ab-initio 
simulation J 1S * 14 * Much attention, in par ticular, has 
been devoted to the characterization of defect^HSEIl i n 
spite of these efforts, the nature and role of defects in dis- 
ordered materials is not fully understood and numerous 
questions remain, especially with regards to relaxation. 
For example, the usual definition of defects — such as va- 
cancies and interstitials — cannot be directly applied to 
disordered or amorphous materials: while 3- and 5-fold 
coordinated atoms are often taken as defect centers in a- 
Si, other medium-range defects, such as string of short or 
long bonds, might also have to be co nsider ed, as shown 
recently by Drabold and collaborator J 32 l 33 [ 

We study here the role of defects with respect to relax- 
ation by examining the response of a-Si to plastic defor- 
mations. In crystalline materials, where the phenomenon 
is well understood, plasticity is attributed to defect nucle- 
ation or dislocation motion^®* in disordered materials, 
the response is harder to define structurally, as demon- 
strated by a number of studies in systems ranging from 
metallic^ and polymeric^ glasses to granular materi- 
als^ and colloidsSS. Helder et a/P^, however, have shown 
that, during irradiation with high-energy heavy ions, a- 
Si deforms plastically in the same manner as conven- 
tional glasses, i.e., defects, irradiation and plasticity are 
directly related in amorphous materials. It is this still 
imperfectly understood relation that motivates the work 
presented here. 

Following initial suggestions by Argon 49 , it has been 
proposed recently that plasticity is caused by collections 
of shear transformation zo nes th at operate as localized 
centers for the deformations) 4 ^ 5 -^^ In a-Si, these plastic 



deformations are attributed to the presence of liquid-like 
particled 52 * 53 * associated to 5-fold coordinated atoms. To 
verify these ideas, elastic and plastic deformations in a- 
Si were investigated by Talati et al. 54 using classical 
MD simulations based on the Stillinger-Weber^ and the 
TersofT 8 potentials; it was concluded that, even though 
the general behavior of the stress-strain curves associ- 
ated with elastic and plastic deformations are similar to 
that for other disordered materials, details as to the na- 
ture of the defect responsible for plasticity depend on the 
particular potential used. 

We revisit this question here using classical MD simu- 
lations and the Environment Dependent Interatomic Po- 
tential (EDIP) 5 6 . Various points defects and their effects 
on elastic constants have been characterized by Allred et 
a/Pausing the EDIP potential 56 . These authors have 
shown that the elastic constants vary in a roughly lin- 
ear fashion with defect concentration up to ~3%, in line 
with experiment, suggesting that EDIP is suitable for in- 
vestigating the elastic and plastic deformations in a-Si. 
Our simulations were performed using the classical MD 
package (LAMMPS^Pl 

Most previous studies of elastic and plastic deforma- 
tions in amorphous materials have been performed at K 
or at low temperature, with a focus on the disordering 
process. Here we examine the situation in a-Si by in- 
vestigating the system's response to variations in shear 
velocity and temperature, as both parameters contribute, 
in their own way, in forcing the system to overcome en- 
ergy barriers and explore the potential energy surface; 
as mentioned previously, struc tura l changes are strongly 
correlated with imposed straiiF^HIl. 

This paper is organized as follows: in the section [ll[ we 
describe the details of our simulation model. In the next 



section (III), we present the results obtained by applying 



mechanical shear deformations at different shear veloci- 
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ties and temperatures. In the section [TVj we present the 
structural analysis as a function of shear velocity. Af- 
ter discussion (section [V]) of our results, we will give our 
main conclusions (section VI). 



II. COMPUTER SIMULATION MODEL 

Classical MD simulations have limitations, but they 
are unavoidable to simulate systems large enough to limit 
the impact of size effects. For the problem at hand, 
charges and electronic effect should not play a direct role 
and therefore a description in terms of suitable classical 
interactions is appropriate. 

The EDIP functional form involves a two-body radial 
term for bond-stretching interactions and a three-body 
angular term for bond-bending interactions, and each of 
these depends strongly on an effective coordination num- 
ber Z. We have chosen this potential for its ability to 
reproduce a wide range of zero-temperature properties 
of Si, including elastic constants, bulk crystal structures 
and point defects 6 ^ — in particular, EDIP describes 
accurately the vacancy formation energy. This potential 
has been used to study t he ion -beam induced amorphiza- 
tion of crystalline silicorl^^ and the crystallization^^ 
of a- Si starting from amorphous-crystalline interfaces 6 - 3 ^. 

The melting temperature of crystalline silicon obtained 
using EDIP is 1500 K and it is 1200 K for a-Si. In both 
cases, this is roughly 200 K below experimental values 
(1685 K and 1420 K, respectively^). Both phases are 
therefore stable or, at least, metastable in the tempera- 
ture range investigated here. 



A. Sample preparation 

Experimentally, a- Si can be obtained by laser- melting 
and quenching, by chemical vapor deposition 6 -^, or by ion- 
irradiation^. Numerically, our models of a- Si are gener- 
ated using the modified Wooten- Winer- Weaire (WWW) 
bond-switching algorithm, which can produce perfect 4- 
fold-coordinated random network d 1 * 67 ! in good structural 
and electronic agreement with experiments 67 . A 1000- 
atom a-Si cell was first constructed, then duplicated in 
order to obtain a cubic box containing 8000 atoms. The 
latter was then annealed in the NPT ensemble with pe- 
riodic boundary conditions over 150 ns at 300 K so as to 
obtain a well relaxed models of amorphous silicon. The 
final structure was found to contain less than 3% of de- 
fects (mostly 5-fold coordinated atoms), the rest being 
perfectly 4-fold coordinated, in agreement with the best 
finite-temperature models available in the literature. 



B. Procedure for mechanical shear deformations 

The mechanical shear deformations on the 8000-atom 
a-Si model are introduced as follows: the cell is first equi- 



librated (at the target temperature — see below) with 
periodic boundary conditions in all three directions, for 
at least 5 ns. Three different regions are then defined, 
as shown in Fig. [T] an upper and a lower "wall" , each 
containing 1000 atoms, and a region of mobile or active 
particles; periodic boundary conditions are now imposed 
only in the lateral (x and z) directions, while surface 
atoms are fixed in the ^-direction. This procedure is sim- 
ilar to that used by Mokshin et a/P^on single-component 
Lennard- Jones amorphous systems. The thickness of the 
walls, about 6 A, is larger than the EDIP cutoff 5 6 , thus 
ensuring that all particles within the active region (40 A- 
thick) share the same physics. 

The mechanical shear deformations are generated by 
moving the walls at fixed shear velocity (shear rate) v s 
along the shear direction (x). The shear velocity (v s ) 
measures the speed of the deformation while the strain 
rate (e) gives the change in strain with respect to time 
and it corresponds to the shear velocity devided by the 
distance between the walls. In practice, the lower-wall 
position is fixed and only the particles in the upper wall 
are moved by imposing a constant displacement v 3 x h at 
every timestep h (1 fs) — the wall particles are otherwise 
frozen in place. The amorphous nature of the walls en- 
sures that crystal growth, if it occurs, is not induced by 
the boundaries as in the case of crystallization studies us- 
ing amorphous-crystalline interface d 69 " 72 !. Note that the 
cell is aged during 1 ns before applying the deformations 
in order to ensure proper relaxation after releasing the 
periodic boundary conditions along y. 

The equations of motion for mobile particles are in- 
tegrated using the velocity Verlet algorithm. Since an 
additional force is imposed in the shear direction (x), 
only the components in the "neutral" y and z direc- 
tions are considered for computing and rescaling the 
velocities to constant temperature during the deforma- 
tion process. The pressure at time t is computed from 
instantaneous atomic positions i = 1,N] and 

forces acting on the particles [^(t), i = 1,N] using 

P(t) = f k B T + 3^ Y^=i ^(*) • with N the num- 

ber of particles, V the volume, T the temperature, and 
fc# Boltzmann's constant. Statistics such as potential 
energy and coordination number were obtained at zero 
shear velocity and will be used as reference for assessing 
the effect of the shear deformations. 

The elastic and plastic deformations of materi als are 
usually analyzed in terms of stress-strain curved 5 -^ 57 * 73 * 
where the strain corresponds to the maximum displace- 
ment in the shear direction with respect to the distance 
between the two walls, which is kept constant. We prefer 
to characterize the deformations in terms of the potential 
energy difference (PED) AE = E p — Eq , where E p and 
Eq are the potential energies of sheared and non-sheared 
systems; Eq is computed before switching on the shear 
deformations. This analysis has the advantage of pro- 
viding direct microscopic information and relates more 
readily to structural changes during shearing. Structural 
properties like radial distribution functions and the co- 
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ordination number will also be used to determine how 
the properties of a- Si are affected by mechanical shear 
deformations. 




FIG. 1: Snapshot of a typical configuration of the 8000- 
particle a-Si-model. Mobile particles are positioned between 
two parallel walls, each containing 1000 atoms, used for ap- 
plying the shear deformations (see text for more details). 



III. SHEAR DEFORMATIONS RESULTS 

In this section, we examine the effects of varying the 
shear velocity and the temperature on the system's struc- 
tural disorder and defect concentration. 



A. Mechanical shear deformations at 300 K 

We fix the temperature at 300 K and vary v s be- 
tween 10 -5 and 8 x 10 -2 A/ps. For all values of v s but 
the smallest, simulations are run until a strain of 20% 
is reached; e.g., for a shear velocity of 10 -2 A/ps, the 
simulation time is 0.8 ns. For the lowest shear velocity 
(v a = 10 -5 A/ps), because of computational limitations, 
we stopped at 12% strain, corresponding to a simulation 
time of 500 ns. In all cases, simulations are long enough 
for the plastic deformation threshold — defined by the 
onset of irreversible deformations — to be reached; this 
occurs in a-Si at a strain of about 10%. The exact value 
of this threshold depends on the shear velocity and we 
found it to increase with v si in agreement with previous 
results (see Ref. 57 , for example). 

It is well known that stress appears immediately follow- 
ing a shear deformation. At short times or small strain 
values, the stress increases almost linearly and saturates 
at high strain values. This constant stress corresponds 



to the sheared steady-state character of plastic deforma- 
tions. This general behavior is common for polymeric 
material s^, metal lic glasses 44 , colloids^ and amorphous 
material J 48 * 52 ^ ^. The response of the system to me- 
chanical shear deformations c an be analyzed, in general, 
in terms of stress-strain curves^ 4 ^. As mentioned above, 
we prefer to use the PED between sheared and non- 
sheared cells, as it provides a direct and simple descrip- 
tion of the microscopic structural deformations that take 
place under shear deformations. 

Figure [2] shows (a) AE, (b) the internal pressure as 
a function of strain for different shear velocities and 
(c) AE higher strain values using a shear velocity of 
8 x 10 -2 A/ps( Arrows indicate the direction of increas- 
ing shear velocity). Evidently, changes in the potential 
energy are a manifestation of structural rearrangements 
that have taken place with respect to the initial, un- 
strained model. Thus, mechanical shear — and therefore 
strain — increases the potential energy, i .e., the disorder, 
of the system (cf. also Ivashchenko et alW^). 

Changes in the potential energy as a function of shear 
take place over two different regimes. For small shears, a 
quadratic behaviour is observed, associated with elastic 
and reversible deformations; the quadratic nature of the 
energy is independent of the rate at which the strain is 
applied, as one can see from the overlap of the curves 
corresponding to different shear velocities. The point at 
which the system crosses over to the second, high-strain 
regime, however, does depend on v s ; the larger v s , the 
longer the elastic regime persists. This behavior is per- 
fectly echoed in the pressure — a larger v s allows the 
system to reach a much higher negative pressure before 
plastic deformations are forced. The potential-energy 
maximum also corresponds to a stress maximum in the 
stress-strain curves or the yield stress (not shown). This 
behavior is common for elastic and plastic deformations 
of polymeric^, metallic glasses 44 , and amorphous mate- 
rials^. 

The onset of plastic deformations is characterized by 
an overall relaxation of the system, as we observe in 
Fig. [2|a). The higher the yield stress (defined by the 
pressure, for example), the larger the relaxation; for the 
smallest v s , 10 -5 A/ps, for example, there is very little 
relaxation after the plastic deformation peak is reached. 
We note however that, after an initial drop, the potential 
energy increases again, but at a relatively slow rate. It 
would eventually saturate to a steady-state level, beyond 
the reach of our simulations, which corresponds to the 
plasticity regime associated with a flowing steady state, 
as one can see in Fig. [2|c). In this simulation, performed 
with a significantly faster shear velocity, the potential 
energy saturates within the run time indicating that the 
system has reached the flowing steady state and therefore 
its maximum disorder. The inset of this figure shows the 
details for strain values up to 1.0. For all shear velocities 
considered here, the mean value of the potential energy 
in the plastic steady state is higher than its value before 
shear deformations are applied. In other terms, the sys- 
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tern exhibits an irreversible plastic deformation leading 
to increased internal strain as in the case of me tallic and 
silica glasses under high energy irradiatio n 174 * 75 !. As for ir- 
radiation, the mechanical shear deformations lead to the 
formation of new defects in a-Si; this will be discussed in 
the section HY Bl 



In the case of plastic deformations, the mean value of 
the steady-state potential energy increases with v s , as 
can be seen in Fig. [2|a) — the curves for different val- 
ues of v s are almost parallel, shifted to lower energies for 
smaller shear velocities. In all cases, we observe a drop 
in the potential energy after the crossover from elastic to 
plastic deformations. This suggests that, at this point, 
the system has to overcome an energy barrier to launch 
a relaxation cascade. The height of this barrier decreases 
with increasing shear, and the probability of crossing it 
at lower strain increases for systems under slower shear 
rates, explaining the observed behavior. Because of the 
Boltzmann factor governing the jump rate, the transition 
is also be pushed to lower strain values with increasing 
temperature. For example, the minimum strain value is 
found to lie between ~2% for the lowest shear velocity 
and ~8% for the highest in our 300 K runs while, accord- 
ing to Talati et a/P^, it can reach 20% at K. Clearly, 
temperature plays an important role on the yield-stress 
value. 
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At 300 K, the PED between sheared and annealed 
a-Si is about 0.02 eV/atom for v s = 10 -5 A/ps and 
0.08 eV/atom for v s = 8 x 10 -2 A/ps: evidently, the 
structural changes taking place in a-Si are strongly cor- 
related to the applied shear velocity. While it was not 
possible to apply a slower strain rate, these results sug- 
gest that the steady-state PED from the annealed state 
would tend to zero as the strain rate decreases. This 
question will be discussed in more details in Section |V[ 



The shear-induced disorder leads to the forma- 
tion of higher-density regions associated with higher- 
coordinated, "liquid-like" particles. This would explain 
the increase in amplitude of the negative pressure ob- 
served in Fig. (2jb) for different shear velocities. For very 
small strain, the pressure decreases harmonically and 
reaches a minimum before settling to a steady-state value 
associated with the plastic-deformation regime, which is 
strain independent (as in the case of the PED). For larger 
strain, the pressure does depend on the shear velocity: 
the change is 2 GPa for v s = 8 x 10 -2 A/ps, while it is 
only 1 GPa for v s = 10 -3 A/ps. The plastic deformations 
are accompanied by a noticeable drop in the pressure am- 
plitude This behaviour does not appear to be caused 
solely by the unusual phase diagram of silicon, where the 
liquid is denser than the solid, as it is also observed in 
metallic glasses and is rather due the concentration of 
defects observed in the plastic regime (see below). 
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FIG. 2: (a) Potential energy difference and (b) pressure dif- 
ference between sheared and non sheared system at 300 K 
as a function of strain, for shear velocities in the range 10 -5 
to 8 x 10 -2 A/ps. The potential energy of the non-sheared 
system (—4.3707 eV/atom) is used as reference to calculate 
the PED. The arrows refer to the direction of increasing shear 
velocities, (c) Potential energy difference for higher strain val- 
ues using a shear velocity of 8 x 10 ~ 2 A/ps. The inset zooms 
on the low- strain values. 
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B. Deformations at fixed shear velocity: 
temperature effects 

In the previous section, we saw that a higher shear 
velocity leads to a larger deformation energy in the plas- 
tic deformation regime at a given temperature. In this 
section, we discuss the role of temperature on the elas- 
tic and plastic properties of a- Si at fixed shear veloc- 
ity Very few such studies have been reported on dis- 
ordered systems and most studies on a- Si were carried 
out at low temperature — or 300 K 76 . The effects of 
temperature and shear velocity in amorphous polymers 
and Lennard- Jones glasses were investigated by Rottler et 
a/P^ who found that they are akin to redefining the time 
scale for structural modifications: these parameters mod- 
ify the properties of the glassy state by altering the aging 
process and inducing rejuvenation. At low temperature, 
atomic diffusion in glassy materials is completely local 
and negligible. Increasing the temperature or the strain 
facilitates the possibility for particles to escape from the 
cage formed by surrounding particles, thus accelerating 
diffusion and relaxation. 

In order to isolate the specific effects of temperature, 
we varied the latter between 10 and 1200 K at fixed 
shear velocity v s = 10 -3 A/ps, during 8 ns for a total 
strain of 20%, allowing the system to age during at least 
5 ns at each temperature. Fig. |3ja) shows the strain 
dependence of the PED as a function of temperature. 
For temperatures up to 900 K, the low-strain harmonic 
regime is again observed, and this is followed by an elas- 
tic to plastic transition leading to a steady state. The 
steady-state PED decreases rapidly with temperature, 
from -0.12 eV/atom at 10 K to -0.04 eV/atom at 300 K, 
to less than 0.01 eV/atom at 900 K. Increasing temper- 
ature is therefore similar to reducing the strain rate. 

The results are different for temperatures above 900 K, 
on which we zoom in in Fig. [3fb): rapid, very small 
(—0.01 eV/atom) energy fluctuations, but no saturation, 
are observed: clear plastic deformations can no longer be 
defined and shear effects are completely compensated by 
thermal relaxation. 

Interestingly, the steady-state PED curves are all 
nearly parallel, a behavior also observed by Rottler et 
a/P^on stress-strain curves at different temperatures. At 
low temperatures, the potential energy increase is essen- 
tially controlled by the shear-imposed displacements and 
we observe almost no self-diffusion; at high temperature, 
in contrast, our results reveal a competition between the 
stress imposed by shear deformations and temperature- 
enhanced self-diffusion that favors annealing. 



IV. STRUCTURAL MODIFICATIONS DURING 
SHEAR DEFORMATIONS 

We discuss in this section the microscopic changes that 
take place under shear, focussing on the radial distribu- 
tion function (RDF), g(r), that provides a global picture 
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FIG. 3: (a) Potential energy difference (PED) as a function of 
strain at fixed shear velocity v s — 10 -3 A/ps for temperatures 
(from top to bottom) 10 K, 100 K, 200 K, 300 K, 400 K, 
500 K, 600 K, 700 K, 800 K, 900 K, 1000 K, 1100 K and 
1200 K. For clarity, the high-temperature curves (800-1200 K) 
are reproduced in (b) on a finer scale. 



of structural changes, and local coordination, which we 
can relate to solid-like or liquid-like behavior. 



A. Radial distribution functions 

Figure [4] (a) shows the average RDF computed for dif- 
ferent strain rates at 300 K. To correct for boundary ef- 
fects, the local, atom-specific RDF is normalized by the 
fraction of the surface of the sphere surrounding it, 4tt r 2 , 
that fits into the system. Thus, the normalization factor 
depends on the distance of the particle from the walls 
in the y direction. Because deformations in the elastic 
regime do not lead to permanent damage, we focus here 
on the plastic steady-state regime and average over the 
last 200 configurations at the maximum strain of 20%. 

Overall, very little difference between the various 
RDFs is observed: the positions of the first and second 
peaks are essentially unchanged and match those of the 
annealed model. While the width of the first- neighbor 
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peak increases slightly, most changes take place near the 
second-neighbor peak (cf. inset): as the height of this 
peak decreases, a new structure appears at shorter dis- 
tances (~2.8 A), which develops as the shear velocity is 
increased. This feature is related to the appearance of 
5-fold coordinated atoms, as already shown in Refs. 52 53 : 
plastic deformations are associated with the presence of 
higher-coordination liquid-like particles, so that the co- 
ordination number increases with shear velocity, which 
confirms the tendency of mechanical shear deformations 
to increase disorder. 

At fixed shear velocity (v s = 10 -3 A/ps), the new 
structure formed by shear deformation is more pro- 
nounced at lower temperatures, as one can see from the 
inset to Fig. [4] (b) where we present the RDFs at differ- 
ent temperatures. As temperature increases, the system 
can more easily compensate for the applied strain, which 
leads to a less pronounced structure between the first and 
second peak, until it vanishes completely at high temper- 
ature, thus confirming that the action of shear dimishes 
as temperature increases. The fact that the amplitude 
of the new structure increases with shear velocity and 
decreases with temperature signals, again, the existence 
of a competition between shear velocity and temperature 
as to their effects on the structure. 



B. Coordination number and defects fraction 

We examine now the evolution of the number of coordi- 
nation defects under shear, counting atoms within a cut- 
off of 2.8 A which corresponds to the minimum between 
the first and second neighbor peaks in the annealed-cell 
RDF. We focus on 5-fold-coordinated atoms since the 
number of other coordination defects remains very small 
during the whole process: the initial system contains 
~97% of perfectly coordinated atoms, less than 3% of 
5-fold and only 0.06% of 3-fold-coordinated atoms. 

The evolution with strain of the population of 5-fold- 
coordinated defects is displayed in Fig. [5] for the various 
shear velocities considered, at 300 K. While the behav- 
ior of these curves is similar to that of the PED and 
the pressure presented earlier, with well-defined elastic 
and plastic regimes, it differs in two notable ways: First, 
there is no decrease in the number of defects at the yield 
stress, but only a sharp flattening of the curve. This 
suggests that, although the concentration of defects def- 
initely increases in the elastic regime, there is no signif- 
icant structural reorganization. At the elastic to plastic 
transition, even though some energetic relaxation takes 
place at the onset of plasticity, the defects created do not 
anneal. Second, the population of defects for the smallest 
shear velocity (10 -5 A/ps) evolves along three regimes: 
a very small elastic region for strains below 0.01, followed 
by a steady increase of the concentration of defects asso- 
ciated with a fall of equilibrium, and finally a transition 
to the steady-state plastic regime at a strain of ~0.03. 
For slow enough shears, therefore, it appears as though, 
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FIG. 4: Radial distribution function (a) as a function of shear 
velocity at 300 K, and (b) as a function of temperature at v s = 
10" 3 A/ps. The arrows indicate the direction of increasing 
shear velocity or temperature (see text). 



in contrast to crystalline systems, the perfectly elastic 
regime disappears: amorphous silicon shows a continu- 
ous distribution of energy barriers leading to structural 
rearrangements, and any amount of shear can move the 
system from one minimum to another^^^. 

Shear affects the structure very significantly. For the 
highest shear velocity considered, 8 x 10 -2 A/ps, the pro- 
portion of 5-fold-coordinated defects increases from as 
little as 3% initially to 23% in the steady state regime; 
it falls to about 12% upon decreasing the strain rate by 
a factor of 80 (v s = 10 -3 A/ps), and is only ~8% for 
the very lowest shear velocity investigated (10 -5 A/ps). 
Thus, even at 300 K, the system manages to anneal it- 
self. Even for the highest defect levels, other defects (e.g., 
3-fold and 6-fold coordinated atoms) remain rare, well 
below 0.5%. Interestingly, these results, obtained using 
the EDIP potential, are in good agreement with those 
obtained by Talati et al. using the Stillinger- Weber po- 
tential. 
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FIG. 5: Fraction of 5-fold coordinated atoms as a function 
of strain. The different curves correspond to shear velocities, 
as indicated. The arrow indicates the direction of increasing 
shear velocity. 



C. Localization of plastic deformations: defect 
distribution 

The propagation of plastic deformations can be stud- 
ied by monitoring the proportion of defects — here 5-fold 
coordinated atoms — throughout the cell. To this end, 
we divide the system into 12 layers along the ^/-direction 
and calculate the concentration of defects as a function 
of strain; the first and last layers (the walls) are not con- 
sidered since they are used for applying the strain. 

The results are shown in Figure [6] for three different 
shear velocities at 300 K. For all cases considered we 
find that, at low strain, in the elastic regime, the distri- 
bution of defects is low but uniform across the system; 
at the elastic-to-plastic transition, however, strong inho- 
mogeneities appear along the ^-direction. At high shear 
velocities, 10 _1 A/ps , Fig.[6](c), 5-fold defects form near 
the center of the system and reach a very high concen- 
tration (>30%) in the steady-state regime: plastic defor- 
mations are highly localized in these central layers and 
the 3 to 4 layers near the walls remain largely unaffected. 
At smaller shear velocities, plastic deformations propa- 
gate through a wider portion of the cell, reaching the full 
width of it for v s = 10~ 5 A/ps [Fig. ^ (a)]. For the lat- 
ter case, at low strain, defects first appear in the layers 
closest to the walls (layers (2 and 11) before propagat- 
ing across the system. In the plastic regime, the defect 
density increases almost linearly with the distance to the 
walls, with all layers affected by plastic deformations. 

A direct comparison of the figures [Figs. [6] (b, c and d)] 
or by computing the distribution of the defects fraction 
across the y-directions, we found that the width of the 
region affected by shear deformations - where the the de- 
fects fraction increases significantly - is almost the same 
for two different systems. According to these observa- 
tions, the localization of the shear deformations in a few 



layers is not size dependent at least for systems larger 
than 8000 particles. 

The effect of temperature can be assessed in Fig. [7] 
where we plot the corresponding results at temperatures 
of 500 K, 700 K and 900 K, using an intermediate shear 
velocity 10 -3 A/ps. While defects are localized in the 
few layers near the middle of the structure at 300 K [Fig. 
[6] (c)], the distribution broadens with temperature: all 
layers are affected by the shear at 700 K and the dis- 
tribution becomes almost flat at 900 K, even at small 
strain. 




FIG. 6: Fraction of 5-fold coordinated atoms at 300 K as a 
function of strain and layer index; (a), (b) and (c) correspond 
to shear velocities (a) 10 -5 A/ps, (b) 10 -3 A/ps, and (c) 
10 _1 A/ps, at 300 K. The configurations are divided into 12 
layers along the ^/-direction; layers 1 and 12 are not shown as 
they correspond to the frozen walls. The figure (d) correspond 
to a system of 16000 particles and a shear velocity of 5 x 
10 -4 A/ps. The number of layers is 22. 




FIG. 7: Fraction of 5-fold coordinated atoms as a function of 
strain and layer index, at (a) 500 K, (b) 700 K, and (c) 900 K, 
for a shear velocity of 10 -3 A/ps. 
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V. DISCUSSION 

In agreement with the results observed on stress-strain 
curves of the e lastic and plast ic deformations of differ- 
ent materialPH 47 | 50 | 52tt 54|57i76] the system's response to 

strain can divided into three regimes. At low strain, the 
system responds with strain-independent elastic defor- 
mations, as observed in the PED and the pressure, for 
example. Although the number of local defects, here 
over-coordinated atoms, increases during this phase, a 
large fraction of these are reversible. The second regime, 
at intermediate strain, depends on the shear velocity and 
is associated with the elastic to plastic transition. Here, 
the absolute values of the PED and the pressure reach 
a maximum, then fall off to smaller values following a 
cascade of bond rearrangements that reduces the strain 
buildup, the drop being larger for faster shear rates and 
lower temperatures. Surprisingly, however, this is not ac- 
companied by a similar drop in the number of defects. As 
the energy is released, the defects are therefore stabilized 
by irreversible atomic repositioning. The third regime 
corresponds to steady-state plastic deformations and oc- 
curs at high strain values. In this regime, the various mi- 
croscopic and thermodynamic quantities are essentially 
strain-independent. For the pressure, the plastic regime 
rapidly reaches a plateau after the break-down. For the 
potential energy, the convergence is somewhat slower and 
the steady-state value is only slightly lower than the max- 
imum at the elastic-to-plastic transition. The curves are 
also shifted to lower energies with decreasing shear ve- 
locity or increasing temperature^. This contrasts with 
the microscopic density of 5-fold defects, which shows a 
steady increase and is not yet converged for most values 
of the shear velocity and temperature at a strain of 0.2. 

Our results on the effects of shear and temperature al- 
low us to better characterize the spatial inhomogeneities 
involved in the response to shear. While defects form 
almost uniformly across the atomic layers in the elas- 
tic phase, plastic deformations are concentrated within a 
narrow region near the center of the box at low tempera- 
tures and high shear rates, and the fraction of liquid-like 
5-fold coordinated atoms in these laye rs reaches almost 
40%, in line with the results of RefsfSEHl Indeed, the cell 
is breaking into two parts with an almost fluid interface 
between them. Such a sharp break is costly however and 
is imposed by kinetic considerations: the relaxation time 
is too short to permit a global response of the system 
to the external perturbation. The system is, in effect, 
frozen, and stress can build up significantly, which can 
only be released by a size- wide rearrangement, creating 
the observed sharp break. This inhomogeneity explains 
why the radial distribution function is only slightly af- 
fected by an average defect concentration of more than 
20%. 

As the temperature rises or the shear velocity de- 
creases, the relaxation time becomes comparable with 
the shear rate and the response to perturbation is spread 
to larger and larger regions, preventing a large build-up 



characterized by a significant increase in potential energy 
or defect concentration. And while 5-fold coordinated 
atoms remain the main type of defects by which this re- 
sponse takes place, it is no longer possible to speak of 
localized centers of deformation, since the whole system 
is affected in the plastic regime. 

The relation between temperature and shear velocity 
can be seen in Fig. [8] which shows the evolution of the 
maximum PED and the mean value of the PED in the 
plastic deformation regime as a function of shear veloc- 
ity and temperature; the data are extracted from Figs. [2] 
and [3} For the PED, we observe almost the same power- 
law relation between the shear- velocity and the energy at 
the yield-stress and in the steady state plastic regime, in 
both cases with an exponent of ~ 0.18 in a wide range of 
shear velocities between 10 -3 A/ps and 8 x 10 -2 A/ps. 
For the smallest shear velocity, the power-law relation 
seems to be less reliable and the effects of shear appear 
to be more important than could be expected. Since ther- 
mal effects are included in the reference potential energy, 
this suggests that there is no critical rate for which ther- 
mal energy completely obliterates the effects of shearing 
even though plasticity will spread to always larger frac- 
tions of the system. 

In the bottom panel of Fig. [8| a similar power-law be- 
havior is observed for temperatures between 150 K and 
900 K, with an exponent of —0.9. Again, even at high 
temperature (but below melting), the effects of shear 
are not fully compensated but are strongly reduced. At 
low temperature, the power-law relation breaks down be- 
cause the bonding energy introduces a hard threshold on 
the amount of strain that can be stored in the system. 

The origin of these power-law relations can be related 
to the energy landscape structure of the system. Because 
the activation energy necessary to cross a barrier is un- 
corrected with the energy difference between the top of 
the barrier and the final minimum, only the shape of the 
forward activation-energy barrier distribution is impor- 
tant^. From numerical calculations 78 , it is known that 
this activation-energy barrier distribution at T = K is 
a continuous function that can be fitted to: 

Gfb(E) = AE exp {-^g^} , (1) 

where E is the barrier energy, A a normalisation fac- 
tor, and E r pB and ctfb two parameters that depend on 
the system and the degree of relaxation. At finite tem- 
perature, for a well-relaxed sample, barriers below ksT 
are already consumed. The shearing here effectively in- 
creases the system's energy, thus decreasing the height of 
available barriers, and, from there, increases the crossing 
probability at a given temperature. Because energy bar- 
riers of any height may exist, it is possible for the system 
to accommodate shearing at least partially, at any tem- 
perature, as long as the shearing rate is slow with respect 
to the attempt-to-jump frequency. 
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using classical MD simulations with the Environment De- 
pendent Inter-atomic Potential (EDIP). The shear defor- 
mations are introduced by moving one wall in the shear 
direction^. We find that both energetic and structural 
properties of the system depend strongly on the shear ve- 
locity and the imposed strain, as well as the temperature 
at which the deformations are applied. At low tempera- 
tures and for all shear velocities investigated, we observe 
a systematic increase in disorder associated with an in- 
crease in the fraction of coordination defects, in agree- 
ment with the results of the reference d 54 * 73 * . Interest- 
ingly, the shear-induced energy can be written in terms 
of a power-law of both temperature and shear velocity: 
AE oc vfT~P, with a ~ 0.18 and j3 ~ 0.9. Finally, we 
observe a very strong dependence of the spatial distribu- 
tion of defects on these two quantities. For low temper- 
atures or high shear velocities, defects are localized in a 
narrow region in the middle of the cell. As the temper- 
ature increases or the shear rate slows down, this region 
becomes wider until it covers the whole system. This ef- 
fect is due to the existence of a continuous distribution of 
activation energy barriers that allows the system to relax 
at any temperature provided that it has ample time to 
do so, as determined by the shear rate. 



FIG. 8: Log- log plot of the potential energy (maximum and 
mean value in the steady state for high strain values) as a 
function of (a) shear velocity at fixed temperature (300 K), 
and (b) temperature at fixed shear velocity (10 -5 A/ps). The 
inset on the panel (b) shows the same curve in a linear scale. 

VI. CONCLUDING REMARKS 

We have investigated the properties of amorphous sil- 
icon subject to external mechanical shear deformations 
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